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A summary of recent calculations of Casimir forces between a collection of N- 
dielectric spheres is presented. This is done by evaluating directly the force on 
a sphere constructed from a stress tensor, rather than an interaction energy. 
A loop integral formulation is also discussed where we rewrite the expressions 
for the force in terms of loop integrals for the effective classical propagation of 
the electric and magnetic fields. 



1. Introduction 

An important step in better understanding Casimir forces and their ap- 
plication to nano-scale environments is the nature and influence of the 
geometries used. If we consider a collection of small bodies and their inter- 
actions, a quantitative understanding is necessary to address any issue of 
applications. One can calculate the forces between bodies directly by using 
a stress tensor evaluated on the body in question. This is an experimentally 
accessible prediction of the theory e.g. three body forces can be measured 
between dielectric spheres 1 for critical Casimir forces. The choice of stress 
tensor in media is not unique and depends on what consistency criteria is 
used. Typically the differences will show up both in the scale of the forces 
and the higher order curvature corrections, and one may view this as a way 
of finding the correct low energy description. 

To address the nature of such interactions between bodies 2-5 Casimir 
interaction energies have been evaluated for collections of compact objects 
interacting with different force carrying fields (electromagnetic and scalar). 
The approach taken has been to evaluate a suitable energy functional in- 
tegral using a T-matrix, whereby an interaction energy can be deduced, 
normalised with respect to their energy when separated at infinity. 

In this talk I summarise recent work we have done on calculating the 
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Casimir force on a single sphere, in an iV-sphere setup. 6 By using a multiple 
scattering approach to evaluate essentially the classical scattering Green's 
function of the configuration, we are able to evaluate the force directly 
on the sphere. The total path length plays a fundamental role and what 
leads one to a loop description. 7 Similar loop descriptions have appeared 
before in 8-10 where a semi-classical type expression of the density of states 
is found, and in 11 a worldline description for scalar fields is given. 

2. The TV-Sphere configuration 

The question we are addressing is how to calculate the force on a particular 
sphere as a result of all the interactions with the remaining spheres in a 
particular static configuration. The Casimir force on a sphere (given by the 
ball B 2 ) in the j-direction due to the effects of the iV-sphere system is given 

by 

F j (l\N-l)= [ d 3 xV l T^{x). (1) 
Jb 2 

The stress tensor we choose is the standard vacuum expression (which is 
consistent with the Lorentz force law 12 ) 

T t3 {x) = ^(x)E J (x)+B l (x)B J (x) i%(|E(x)| 2 + \B(x)\ 2 ), (2) 

where x G B 2 and it is understood that we are taking the limit for the 
initial and final points. We then need to evaluate the scattering correlation 
functions (whilst dropping the direct modes of propagation) 

limE 1 (x)E J ( 1 /)= f f dudu>'<p? a {x-,uj)*Ei t {y-,u>% (3) 

y^ x Jo Jo 

and similarly for magnetic fields. To construct the scattering two point 
function we write the fields in a mode decomposition 13 of spherical vector 
wave functions that are centred on each sphere centre. Then by applying 
the standard continuity equations at each of the spheres surfaces, one can 
calculate the out modes in terms of the in modes and scattering (Mie) 
coefficients. Assuming that the background in which we are evaluating this 
is filled with quantum noise such that the noise-current two point function 
is non-zero we find for the ./V-body force on a sphere (suppressing the SO(3) 



November 9, 2009 9:9 



WSPC - Proceedings Trim Size: 9in x 6in BabingtonQFEXT09vl 



3 



/ dujkcofh{h^/K B T){l\[a l {wR[l]) 
4?r Jo 

N N 

x ]T ^ M (r[M) • a'iuRii]) ]T A^(r[i,j}) ■ oP 

i=2 j=2 
N 

x^V rb - 1] ^ 1 (r[i,l])]j(fci?[l])/i+(fci?[l])W(a; J R[l])|l). 
-(-D^E^tnf^^ (^) z i<x,AW]. (4) 

Here, a l (wii[i]) are the Mie scattering coefBcicnts in the SO (3) basis for 
sphere i with radius A 1 ^ (r[i, j]) are the translation matrices mapping 
the TE and TM vector wave functions between spheres i and j; the vec- 
tors |1) give the truncation in the L angular momentum quantum number 
(leading to a multipole type expansion); and the W is just the collection 
of the four different contributions that make up the stress tensor, together 
with the two necessary Bessel functions evaluated on the spheres surface. In 
the last line a Wick rotation to imaginary frequencies has been performed 
and the Z-function has been defined for later reference. 

Note the explicit form of the translation matrices involve exponentials 
of the inter-sphere separations 6 and thus it is the total path length that 
plays the key role in understanding the variables of the system. For simple 
setups (e.g. two and three sphere systems) we can evaluate Equation (4) 
in different perturbative regimes e.g. retarded or non-retarded limits using 
static values for the permittivities (see 6 for explicit evaluations and force 
plots of two and three sphere systems). 

3. Loop integrals 

Following on from the observation that Equation (4) involves the total path 
length in the form of a loop, in 7 an attempt is made to develop this further. 
It is similar to the path integrals used in 8 and 11 where a fictitious time 
is introduced as well as mass scale set equal to unity and an appearance 
of Planck's constant in the particle action. As pointed out in, 10 it is mis- 
leading to call this a semi-classical evaluation because of the absence of the 
dimensionful Planck's constant. Concurrent with this is the absence of a 
mass or length scale with which to define a dimensionless action. 



indices) 

F[l|iV- 1] = 
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In 7 a symbolic expression for the Z-function was found that features 
the loop structure in an explicit manner 



zM[Q,x] = ^(F^|exp [- I dqW^n) 



exp 



-s R (c x ,r,n) 



{IJ,FJ} 



(5) 



where x is the initial and final point of the loop, and |I) and |F) are the ini- 
tial and final states (i.e. boundary conditions imposed on the eigenfunctions 
used to represent in and out modes). The generator £>i of translations and 
the loop C x implicitly depend on the background potentials. The connection 
r is formed from the background potentials in which the field propagates. 
One now needs to find a representation of this object. In fact it can be 
given a path integral representation, albeit a classical one i.e. no h fea- 
turing anywhere. One first needs to invert the Hclmholtz operator (here 
partial derivatives have been promoted to covariant derivatives w.r.t. the 
permittivity and permeability) 



A(z,y) = (z\ [- C 2 /Sl 2 V f AV,A-l] : 



\y) 
\y)- 



(6) 



Introducing a world line metric e to implement the Hclmholtz equation on 
physical states, and an integration over the V, operators (again acting on 
physical states), together with an integration over paths one finds for the 
scalar version of the Hclmholtz operator (i.e. two potentials but no spatial 
indices) 

z[n,x}= firV(F x |r/ [^VtopV^Ml 

Jo n=1 J q [0]= q [l]=x 



[dV][de] exp -n / dt[q ■ V + e(V e • V M + 1)] |I X ), 



(7) 



and 



f 

Jo 



dnz[n,x] - z[n,x\. (8) 

loops 

If we perform first the integral over the world line metric, the Helmholtz 
equation is implemented. Performing the integral over the derivatives re- 
turns a configuration space path integral which would require gauge fixing 
the world-line metric. The classical equations of motion then lead to closed 
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geodesies defining the loops. Integration over the frequency provides part 
of the sum over paths, whilst the tau integration gives the windings of the 
loops. 

4. Conclusions 

In this talk I have summarised recent work we have done on calculating 
Casimir forces between spheres using a multiple scattering approach. The 
total closed path length plays a key role in understanding the calculated 
forces and leads to a loop description. By considering the origin of the 
translation coefficients, together with a path integral representation of the 
Hclmholtz operator, one is able to reformulate the loop integral as a sum 
over of all possible loops. 
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